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ABSTRACT 

Very long baseline interferometry (VLBI) observations of SiO masers in Orion Source 
I has enabled for the first time to resolve the outflow from a high-mass protostar in 
the launch and collimation region. Therefore, Source I provides a unique laboratory to 
study mass-loss and mass-accretion in a high-mass protostar. We numerically simulate 
the dynamics of the disk-wind inside 100 AU from Source I. This enables us to inves- 
tigate the balance of different forces (gravitational, magnetic, thermal) regulating gas 
dynamics in massive star formation. In this work, we adopt magnetohydrodynamic 
(MHD) disk-wind models to explain the observed properties of the disk-wind from 
Orion Source I. The central source is assumed to be a binary composed of two 10 Mq 
stars in a circular orbit with an orbital separation of 7 AU. High resolution ideal MHD 
wind launching simulations (which prescribe disk as a boundary) are performed using 
the PLUTO code. The simulations are allowed to run until a steady state is obtained. 
MHD driven disk- wind provides a consistent model for the wide-angle flow from Source 
I probed by SiO masers, reproducing the bipolar morphology, the velocity amplitude 
and rotational profile, the physical conditions, and the magnetic field strength. 

Key words: (magnetohydrodynamics) MHD - methods: numerical - ISM: jets and 
outflows 



1 INTRODUCTION 

In the formation of low-mass stars, disks and jets play a ma- 
jor role in controlling mass accretion and removing excess 
angular momentum. Specifically, low-mass protostellar jets 
are expected to be launched and collimated via magneto- 
centrifugal processes at distances <<10 AU from Young 
Stellar Objects or YSOs (e. g. Konigl & Pudritz 2000). How- 
ever, the balance of forces in proximity to high-mass YSOs 
has been difficult to infer because high extinctions, cluster- 
ing, and large distances hinder attempts to resolve the gas 
structure and dynamics at small radii (<100 AU) where out- 
flows are launched and collimated from accretion disks. 

Radio source I in Orion BN/KL is the nearest exam- 
ple of a high-mass YSO. The close distance and exceptional 
richness in molecular maser tracers enabled the best detailed 
picture of gas dynamics in a massive YSO on scales compara- 
ble with the Solar System. In particular, detailed imaging of 
vibrationally-excited states SiO masers with the Very Long 
Baseline Array (VLBA) enabled resolving the surface of an 
edge-on disk as well as a wide-angle outflow in the collima- 
tion and launch region at radii 10-100 AU (Matthews et al. 
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2010; Figure 1). Additionally, at larger radii (100-1000 AU), 
the wide-angle wind appears to collimate in a straight bipo- 
lar flow, as traced by H 2 and z/=0 SiO maser and thermal 
emission (Plambeck et al. 2009; Greenhill et al. 2012). Mea- 
surement of the 3D velocity fields of SiO masers provides 
a clear evidence of rotation in both the disk and the wind 
along the same axis. Cunningham et al. (2005). explored the 
importance of stellar winds from a luminous central source 
and proposed the idea that SiO masers in Source I may 
arise from the interaction of a spherical wind with an in- 
falling envelope. However, Matthews et al. (2010) argued 
that such a model could not explain in detail the gas dy- 
namics and physical conditions as revealed by SiO masers 
and alternately proposed that the SiO masers arise from 
material ablated from the surface of an accretion disk and 
discussed qualitatively the applicability of various classes of 
disk- wind models to Source I. They find evidence to sup- 
port that magnetic fields may play a role in driving and/or 
shaping the disk-wind from Source I, based on details of gas 
dynamics on small scales, like curved and helical trajecto- 
ries traced by maser motions. Assuming a magnetic field is 
present, a magneto-centrifugally launched wind (e.g., Bland- 
ford k Payne 1982; Konigl k Pudritz 2000; Fendt 2006) 
would be a natural candidate for powering a disk- wind from 
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Figure 1. Velocity field (in km/s) of the SiO v=l,2 J=l-0 maser 
emission in Orion Source I as observed with the VLBA (Matthews 
et al. 2010). The stellar velocity is ~6 km s _1 . The bulk of the SiO 
emission is located within four "arms" of an X-shaped pattern, 
while a ~14 AU thick band with no SiO emission harbors the 
7mm continuum source (Reid et al. 2007; Goddi et al. 2011). The 
continuum emission traces an ionized disk while the SiO masers 
probe material ablated from the surface of the disk in a wind. 




Figure 2. Cartoon Model for Source I. 



Source I. The latter thus becomes an ideal laboratory to 
extend the MHD models from solar-like stars to massive 
YSOs. Recently, Vaidya et al. (2011) explored the interplay 
of radiation forces and magnetic fields to show that radia- 
tive line driving from luminous stars with masses > 30 Mq 
does modify the dynamics of outflows. However, since the 
luminosity for a 10 Mq binary Orion Source I would be less 
than a single 20 Mq , we adopt that model in the ideal MHD 
regime (neglecting radiative force) to numerically simulate 
the dynamics of the disk- wind inside 100 AU. 



2 MODEL ASSUMPTIONS 

The initial conditions adopted to model the disk-wind from 
Orion Source I are reported in Table 1 and are derived from 
either observations and/or modeling of the system. Source 
I is believed to be an eccentric, hard binary with a total 
mass of 2 x 10 M , based on studies of the dynamics of 



YSOs and molecular gas in Orion BN/KL, as well as N- 
body simulations (Goddi et al. 2011; Bally et al. 2011). 
In particular, proper motion measurements of radio sources 
provided strong evidence that a dynamical interaction oc- 
curred about 500 years ago between massive YSOs Source I 
and BN (e. g., Goddi et al. 2011). N-body numerical simu- 
lations showed that the dynamical interaction between a bi- 
nary of two 10 Mq stars (Source I) and a single 10 Mq star 
(BN) may lead to simultaneous ejection of both stars and the 
BN/KL outflow as well as binary hardening, while preserv- 
ing the original circumbinary disk around Source I (Bally 
et al. 2011; Moeckel & Goddi 2012). Mechanical energy con- 
servation (and virial theorem) implies an orbit semi-major 
axis do = 2 — 5 AU, depending on the estimate of the ki- 
netic energy of the runaway stars and the molecular gas 
flow (Goddi et al. 2011; Bally et al. 2011). For simplicity, we 
neglect the eccentricity and assume that the system is an 
equal-mass circular binary, comprised of two 10 Mq stars 
orbiting in a circle of radius ao = 3.5 AU. 

A high-resolution study of Ha and [O IJA6300 line pro- 
files in a pre-main-sequence spectroscopic binary has shown 
evidence of a bipolar jet as being launched by the whole 
binary as opposed to each star individually (Mundt et al. 
2010). Similarly, we assume that the binary members orbit 
inside the inner radius of a truncated circumbinary disk and 
that the wind is launched at this radius. The cartoon figure 
representing our model assumption is shown in Fig. 2. 

To study the wind launching from Source I, we perform 
axi-symmetric ideal MHD simulations in cylindrical coordi- 
nates (r,0,z) using the PLUTO code (Mignone et al. 2007). 
We essentially prescribe a hydrostatic corona (density and 
pressure) threaded with force- free magnetic field. An adia- 
batic equation of state, P — ^^p 7 , with 7 = 5/3, relates the 
pressure with the density in the flow. For the central grav- 
ity, we treat the system of equal mass binary members as an 
equivalent 5 Mq (reduced mass) object at their bary center. 
The accretion disk which is treated as a boundary forms the 
launching base of the wind (e.g., Ouyed & Pudritz 1997) 
Further, we choose the wind launching point to be at a dis- 
tance equal to the radius of the binary orbit, i.e. ro = ao, 
where the material is rotating with a sub-Keplerian veloc- 
ity Vcj) = x v o- Further, the disk is considered to be denser 
than the hydrostatic corona and this density contrast is pre- 
scribed by 5. In a realistic (thick and hot) disk, the central 
gravity is balanced radially by contributions from thermal 
pressure and centrifugal rotation. This can be ensured by 
consistently deriving the sub-Keplerianity(x) based on the 
choice of density contrast. For example, smaller values of 
5 implies a hotter (denser) disk and a larger contribution 
of thermal pressure resulting in the disk to rotate with sub- 
Keplerian speeds to maintain radial balance. We chose 5 ~ 3 
as a reference value implying the underlying disk to rotate 
at a speed 0.8 times the Keplerian speed. 

The magnetic field strength is prescribed by the value 
of plasma beta (the ratio of thermal pressure to the mag- 
netic pressure) at ro, /3q. In order to put constraints on the 
magnetic field strengths in Source I, we study the wind prop- 
erties with varying f3o (see Sect. 3.3). Finally, we assume the 
number density at ro to be no = 3 x 10 10 cm -3 . Reid et al. 
(2007) argued the disk in Source I to be ionized with a typi- 
cal electron density of ~ 10 7 cm -3 (their "Model C"), which 
corresponds to an ionization fraction of ~ 10 -3 . 
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Table 1. Parameters chosen for the reference run to study MHD 
driven outflow from Orion Source I. 



Quantity 


Reference Value 


Masses of Binary [mi, 7712] 


1OM ,1OM 


Orbital Separation [2ao] 


7AU 


Reduced Stellar Mass [/x*] 


5M 


Wind Launching radius [ro] 


3.5 AU 


Kepler ian Velocity at ro [vo] 


35.7 km s- 1 


Volume Density at ro [po] 


5.0 x 1CT 14 gem" 3 


Plasma Beta at ro \fto] 


1.0 


Density Contrast [S] 


2.85 



The flow is freely allowed to leave the outer boundaries 
(r,z = 50 ao, 150 ao) and symmetric boundary conditions 
are applied on the r=0 axis. The poloidal component of the 
magnetic field at the outer boundaries is set to avoid ar- 
tificial collimation by ensuring zero current (V x B = 0; 
Ustyugova et al. 1999; Porth & Fendt 2010). The critical 
boundary at z=0 mimics the inflow of gas and magnetic 
flux from the disk surface. Here, we ensure the disk to be 
in pressure equilibrium with the medium but with a higher 
density (i.e. colder than the medium). The angular speed of 
the magnetic field line in the disk is equal to the Keplerian 
value and the poloidal velocity is updated each time with 
the value at the grid cell just above the boundary. Further, 
ideal MHD approximation permits to set the flow velocity 
along the poloidal magnetic field (v p \\ B p ). Thus, the mass 
flux at the disk boundary is not fixed but consistently de- 
rived by causal interaction between the outflowing gas and 
the boundary value. Detailed description of these boundary 
conditions are given in Vaidya et al. (2011). 



3 RESULTS AND DISCUSSION 
3.1 MHD driven wind profiles 

At the onset of our reference run, material at the base of the 
disk is ejected in vertical direction along the magnetic field. 
This material propagates outwards producing bow shocks in 
the ambient medium. As this shock propagates outside of the 
domain, the flow collimates due to magnetic forces until a 
steady equilibrium is attained, beyond which no further evo- 
lution of magnetic field lines takes place. The steady-state 
values of the density, poloidal and azimuthal velocities of the 
wind are shown in Figure 3, along with the magnetic field 
lines (in red), the poloidal velocity versors (white arrows), 
temperature contours and MHD critical surfaces. The mag- 
netic field lines ideally mark the morphology of the outflow. 
The left panel of Figure 3 shows the steady-state number 
density of the fluid, whose values range from 10 8 cm -3 to 
10 n cm -3 . Overlaid on the density map are the tempera- 
ture contours, with values ranging from 5000 K closer to the 
central source to 100 K at large distances. The high den- 
sity and temperatures close to the axis reflect the fact that 
the initial hydrostatic conditions are fixed at the boundary 
inside of the wind launching radius (this ensures numeri- 



cal stability of the flow very close to the axis). The middle 
panel of Figure 3 shows the poloidal wind velocity, whose 
values range from 8.0 to 32.0 km s _1 . The lines over-plotted 
show three critical surfaces in a steady-state MHD flow, viz. 
the magneto-sonic surface (yellow), the magneto-fast sur- 
face (black), and the Alfven surface (red). The flow starts 
sub-magneto-sonically from the disk and is accelerated to 
super fast speeds mainly via Lorentz forces. The Alfven sur- 
face marks the locus of points in the flow where ideally the 
kinetic energy of the flow is converted to magnetic energy. 
This results in bending of the field lines thereby generating 
a toroidal field and corresponding hoop stress which colli- 
mates the flow. The right panel of Figure 3 finally shows the 
azimuthal component of the wind velocity, with values rang- 
ing from 10 kms -1 at the outer edge to 25 kms -1 close to 
the inner wind launching point. The transfer of angular mo- 
mentum from the wind to the medium may impart rotation 
to the latter, as observed in SiO masers (Sec. 3.2). 

3.2 Comparison between MHD wind models and 
wind properties observed from SiO masers 

We propose that the SiO masers in Orion Source I are 
excited as a moderately rotating disk-wind interacts with 
the ambient medium via slow MHD shocks and transfer 
of angular momentum. Similar mechanisms have been suc- 
cessfully applied to explain the excitation of H2O masers 
(Kaufman & Neufeld 1996). Here, we qualitatively compare 
the steady- state properties of a single-fluid, rotating MHD 
wind (launched from an underlying accretion disk) with the 
observed physical and dynamical properties of the outflow 
molecular component probed by SiO masers in Source I. 

Three physical properties of the wind from the MHD 
modeling are consistent with SiO maser observations: tem- 
perature/density profiles, morphology, and velocity field. 

The conditions for optimum excitation of the SiO maser 
emission are n(H2)^ 10 9 cm -3 and T< 2000 K, respectively, 
as derived from a radiative transfer analysis by Goddi et al. 
(2009). In our simulations, the regions close to the in- 
ner launching point (z,r < 15 AU) have values of densities 
(10 10-11 cm" 3 ) and temperatures (2000-5000 K) beyond the 
critical values required for SiO maser excitation (Fig. 3, left 
panel). This is consistent with the 14 AU thick "dark band" 
seen in Fig. 1 where no masers are observed. The interme- 
diate region (15 AU < z,r < 60 AU) has typical densities 
of 10 9 cm" 3 and T ~ 1000 - 2000 K (Fig. 3, left panel). 
Here, the wind can provide a strong IR background source 
for pumping the SiO masers in the shocks produced by its 
interaction with the ambient molecular medium. 

In the steady-state, the magnetic collimation and ac- 
celeration result in an open funnel shaped structure of field 
lines (Fig. 3). Such a configuration of magnetic flux sur- 
faces naturally results in a disk- wind with wide-angle bipolar 
morphology, which resembles the X-shaped wind structure 
around Source I revealed by SiO maser emission (Fig. 1). 

Regarding the velocity field, the dynamical properties of 
molecular outflows can be extrapolated from the single-fluid 
MHD wind by analyzing the momentum transfer (Masson & 
Chernin 1993). In this study, the density of the MHD driven 
wind exceeds the mean density of its surrounding medium. 
Such dense and cold winds carry a large amount of momen- 
tum and impart a significant fraction of it to the ambient 
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Figure 3. Number density (left), poloidal velocity (middle) and azimuthal velocity (right) of the reference run in steady-state. Temper- 
ature contours (black) for T = 500 K, 1000 K, 2000 K and 5000 K are overlaid in the left panel. The middle panel shows the positions of 
three MHD critical surfaces - magneto-sonic surface (yellow), Alfven surface (magenta), magneto-fast surface (black). White arrows in 
the right panel represents the velocity vectors which are parallel to magnetic field lines shown by solid red lines. 
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Figure 4. Radial profiles of the magnetic field [in mG] at z = 
20 AU [solid] and z = 60 AU [dashed] for the reference run, ob- 
tained after the MHD wind has reached a steady-state. 

medium, resulting in the shocked material to have veloci- 
ties of the same order as the wind itself (Masson & Chernin 
1993), i.e. in the range of 10 - 30 krns" 1 (Fig. 3, middle 
panel). Consistently, Matthews et al. (2010) measured the 
3D velocities of the molecular gas flows around Source I, to 
vary in the range from 5 km s _1 to 25 km s _1 , with an av- 
erage of 14 km s _1 . In addition, SiO masers show rotation 
seen as a velocity gradient in the radial velocities of the four 
arms of the "X" tracing the disk-wind, which is consistent 
with the trend obtained for V<f, (Figure 3, right panel). 

Therefore, the simulated MHD winds can drive molecu- 
lar flows by interacting via shocks with the ambient medium. 

3.3 Magnetic field strength 

In our models, the magnetic field strength is governed by 
/3q. We have carried out simulation runs with three differ- 



ent values of /3q viz. 1.0, 3.0 and 10.0. Figure 4 shows the 
radial profiles of magnetic field strength at z = 20 AU and 
z = 60 AU for different values of f3o and 5= 2.85. We find 
that the magnetic field strength varies in the range 0.2-0.8 
G (depending on flo) and decreases with radius and/or ver- 
tical height. From this study, we observe that the changes 
in density and temperature are not significant with vary- 
ing f3o (e.g., AT ~ 100 — 200 K when /3q increases by ten 
fold). On the other hand, poloidal and azimuthal velocities 
decrease substantially with increasing values of /3q. Since 
lower field strengths imply lower Lorentz force, the flow is 
not efficiently accelerated resulting in lower velocities. The 
typical velocities obtained from f3o > 3.0 are not enough to 
explain the observed velocities in Source I. Therefore, we as- 
sume as a reference the magnetic field strength correspond- 
ing to 0o = 1, -0.35 G (at r, z =40 AU, 20 AU). Assuming 
equipartition, Matthews et al. (2010) quote a magnetic field 
strength of — 0.3G for Source I, based on kinematics of SiO 
masers. This value is consistent with a measurement of OH 
1665 MHz masers across a — 10 4 AU region around Source 
I, which provided field strengths of 1.8 to 16.3 mG (Cohen 
et al. 2006). Since the OH masers are believed to arise from 
material with n(H2) ~ 10 7 cm -3 (Gray et al. 1992) and the 
magnetic field strength is expected to scale as the square 
root of the gas density, a simple extrapolation from the large 
scales (probed by OH masers) to the small scales (probed by 
SiO masers) gives field strengths of order the value predicted 
by the equipartition assumption. These empirical estimates 
lie comfortably within the range predicted by our model, but 
direct measurements via the Zeeman effect could provide a 
more reliable comparison between models and observations. 

3.4 Model Limitations 

The main qualitative features of a steady MHD driven disk- 
wind are consistent with the physical and dynamical prop- 
erties of the molecular flow probed by SiO masers. A more 
quantitative analysis requires however further theoretical 
and observational studies. In this section, we list the lim- 
itations of the present model and the scope for future work, 
(i) The assumption of treating the disk as a bound- 
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ary allows us to follow wind dynamics till a steady state is 
reached, but neglects the impact of disk dynamics on wind 
propagation. A more realistic model should allow both disk 
and outflow to evolve simultaneously (e.g., Casse Sz Kep- 
pens 2002; Sheikhnezami et al. 2012). Also, we neglect any 
feedback on the disk- wind from a stellar wind (Source I has 
however a rather modest Lboi ~ 10 4 L©; Gezari et al. 1998). 

(ii) In our ideal MHD model, we have assumed a single 
fluid MHD driven wind and have ignored effects due to diffu- 
sion of neutral molecules across the magnetic field. Though 
the regions close to the star are hot and potentially have a 
high degree of ionization, consistent with an ideal approx- 
imation, a more accurate model should incorporate effects 
due to diffusion between the neutral and ionized medium. 
Future simulations will relax the ideal MHD approximation 
and explore effects of magnetic diffusivity. 

(iii) We have ignored any feedback due to molecular 
chemistry. Ideally, a direct comparison between MHD mod- 
els and maser observations would require a detailed multi- 
fluid chemical model. Recently, conditions for molecular sur- 
vival in MHD disk- winds were computed using chemical evo- 
lution of a self-similar MHD wind model (Panoglou et al. 
2012). However, such a chemical model with micro-physics 
of maser excitation is quite complex. As a follow-up, we plan 
to include non-equilibrium cooling and separately follow the 
evolution of neutral and ionized gas fraction in the flow. 

(iv) Including in the modeling the elements discussed 
above, requires measurement of relevant physical quantities 
within the launching region (< 100 AU), i.e. ionization de- 
gree, molecular abundances, and magnetic field strengths. 
High resolution and sensitive observations of Hydrogen ra- 
dio recombination lines (with JVLA and ALMA) will help us 
constrain the dynamical and physical properties of ionized 
gas. Further, radiative transfer analysis of rotational molec- 
ular lines (from ALMA) can constrain molecular gas prop- 
erties and abundances. Finally, polarimetric observations of 
SiO masers (with VLBI) will enable to measure magnetic 
field strength and orientation at the disk/jet interface. 



4 SUMMARY 

We have explored the plausibility of MHD driven wind 
model to explain the disk-wind properties in Orion Source I 
as probed by SiO masers. In particular, we have performed 
axisymmetric MHD wind launching simulations treating 
Source I as an equal mass binary and choosing initial pa- 
rameters derived from observations. We propose that the 
SiO masers are excited as the MHD driven wind interacts 
with the ambient molecular medium in form of shocks. Four 
elements of consistency support this hypothesis: 

(i) Our model consistently reproduces the bipolar mor- 
phology of the disk- wind traced by the SiO masers. 

(ii) The poloidal and azimuthal velocities from the simu- 
lations are in the range of 10 - 30 km s~\ which agree very 
well with the observed average 3D velocity of 14 km s _1 and 
with the rotational component of SiO masers. 

(iii) The physical conditions in the modeled wind are con- 
sistent with the requirements for pumping SiO masers. 

(iv) We predict field strengths of 0.2-0.8 G within 100 AU, 
in agreement with empirical estimates from SiO masers. 



We conclude that the main dynamical features of the 
SiO masers disk-wind in Orion Source I can be explained by 
the MHD wind model. In a broader context, Source I has 
provided a unique laboratory to extend the standard MHD 
models from low-mass to high-mass YSOs, indicating that 
in both regimes similar physics may describe the transfer of 
matter and angular momentum between disks and outflows. 
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